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Abstract. Using a population number synthesis code with detailed binary evolution, we calcu- 
late the distribution of the number of type la supernovae as a function of time after starburst. 
This is done for both main progenitor scenarios (single degenerate and double degenerate), but 
also with various evolutionary assumptions (such as mass transfer efficiency, angular momentum 
loss, and common envelope description). The comparison of these theoretically predicted delay 
time distributions with observations in elliptical galaxies then allows to constrain the evolution- 
ary scenarios and parameters. From the morphological shape of the distributions, we conclude 
that all supernovae la cannot be produced through the single degenerate scenario alone, with 
the best match being obtained when both scenarios contribute. Within the double degenerate 
scenario, most systems go through a phase of quasi-conservative, stable Roche lobe overflow. 
We propose stellar rotation as a possible solution for the underestimation of the observed abso- 
lute number of events, as is the case in many theoretical population synthesis studies. A brief 
comparison with these other studies is made, showing good correspondence under the nontrivial 
condition of equivalent assumptions. We also investigate the influence of different supernova la 
progenitors and evolutionary parameters on the theoretical distribution of the iron abundance 
of G-type dwarfs in the Galactic disk. These stars are good indicators of the entire chemical 
history of the Galaxy, and their predicted metallicity distribution can also be compared to the 
observational ones. This again limits the number of acceptable combinations of assumptions. 
Supporting previous results, the best correspondence is found in the case where both the single 
and double degenerate scenario contribute. 
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1. Introduction 

To distinguish between the single (SD; accretion by a white dwarf (WD) from a late 
main sequence (MS) or red giant (RG) companion) and double (DD; the spiral-in and 
merger of two WDs due to gravitational wave radiation emission) degenerate scenario 
for type la supernova (SN la) formation, we use an updated version of the population 
number synthesis (PNS) code by De Donder & Vanbeveren (2004). This code contains 
detailed binary evolution, without the use of analytical formalisms. SD progenitors are 
assumed to be as given by Hachisu et al. (2008), including the mass stripping effect with 
its strength parameter C\ which can lie between (turning the effect off) and 10. For the 
DD scenario, it is either assumed that every WD merger exceeding 1.4 Mq results in such 
an event, or that it is additionally required to be a merger of two C-0 WDs. Obviously, 
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these assumptions result in different absolute SN la rates, however, the observed trends 
and thus the conclusions are not affected. A parameter study is undertaken, first and 
foremost concerning the fraction j3 of Roche lobe overflow (RLOF) material that is being 
accepted by the accretor star. If f3 < 1, indicating non-conservative RLOF, it is assumed 
that the lost matter leaves the system through a circumbinary disk, therefore taking with 
it the specific angular momentum of L^. If a common envelope (CE) phase is encountered, 
it is treated with the ev-formalism by Webbink (1984). 

There are two typical channels that can lead to a DD SN la: a system can either evolve 
through a stable RLOF followed by a CE phase (typically resulting in a delay time of a 
few Gyr) or through two successive CE phases (with a typical delay of only a few tens or 
hundreds of Myr) . It is important to remark however that the former channel only works 
if the RLOF phase is assumed to be (quasi-)conservative, otherwise systems will already 
merge during this phase and thus not result in a SN la. As in most other PNS studies, we 
find an absolute SN la rate that lies too low compared to observations, in our case by a 
factor of 2-3 at 11 Gyr. While this may be partially due to errors in the conversion factor 
between observational (SNuK) and theoretical (SNuM) units (see also the contribution 
by Pritchet), we present a suggestion for a more physical solution to this problem. There 
are indications that stars born in binaries rotate faster on the ZAMS than normal single 
stars. If this is so, they will also have heavier MS convective cores, which will (even if 
the stars are synchronized later on) result in heavier WD remnants and thus more WD 
mergers exceeding 1.4 M . 



2. Delay time distribution 

Figure 1 shows the theoretically predicted SN la delay time distribution (DTD) for 
different scenarios and values of /3, as well as two observed DTDs. It is obvious that the 
SD DTD (shown for c\ — 3) drops away too fast and too soon in order to keep matching 
the observations at later times. The DD DTD is able to match the observations, but only 
in the case of (3 — 1. This implies not only that (3 indeed needs to lie close to 1, but also 
that those systems which no longer result in a SN la if f3 = (about 80% of the total) 
have gone through such a quasi-conservative RLOF phase, and not two successive CE 
phases. Most of the SD events issue from the WD+MS channel, not WD+RG. The best 
match between the observed and predicted DTD is obtained when we assume that both 




Figure 1. DD DTDs for /3 = 1 (solid black) and j3 = (dashed black), as well as SD DTD 
(dotted gray). Observations by Totani et al. (2008) (black circles) and Maoz et al. (2011) (gray 
squares) . 
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the SD and DD channel contribute, the former with a slightly less high value (ci = 1) 
for the mass stripping parameter than is standard. 

A note on the differences in models obtained by different groups performing PNS 
studies (see also the contributions by Toonen, Nelemans and Claeys) . A comparison study 
is going on to look for the causes of differences in the predictions. It should finally result in 
an explanatory paper. The most important discrepancies are being caused by differences 
in assumptions, mainly about mass and angular momentum loss and the treatment of CE 
phases. Homogenizing the assumptions as much as the codes technically allow results in 
converging models. Some (more minor) discrepancies remain. This seems inevitable given 
the different treatment of the "single star tracks" , i.e. the stellar evolution of individual 
binary components. For example, the use of the Hurley et al. (2002) prescription, versus 
using full binary evolution. 



3. G-dwarf metallicity distribution 

A second distribution which allows to constrain SN la progenitors is the G-dwarf metal- 
licity distribution. Because of their very long lifetime, G-type dwarfs in the Galactic disk 
arc excellent indicators of the entire chemical history of this region. The metallicity 
([Fe/H]) distribution of these stars is critically affected by the SN la rate throughout the 
Galaxy's history, and thus also by the progenitor assumptions. The thus predicted distri- 
butions can be compared to observed distributions for a cylindrical solar neighborhood, 
in this case those by Holmberg et al. (2007). We use an updated version of the chemi- 
cal evolution model used by De Dondcr & Vanbeveren (2004) for a similar study. Other 
studies have been performed in the meantime by e.g. Greggio et al. (2008) and Matteucci 
et al. (2009) (both not with a full PNS model but with various adopted DTDs), as well 
as very recently by Kobayashi & Nakasato (2011) (for SD progenitors only). We use a 
binary fraction (the fraction of stars with a stellar companion) of 70%, which is required 
to be that large in order to yield sufficiently high SN la rates. For the galaxy formation, 
we assume either the two-infall model or a constant star formation rate (SFR). 

In case of the SFR obtained from the two-infall model, and for contribution by both 
the SD (with c\ = I) and DD channel, we find a metallicity distribution which peaks at 
approximately the right location, however it is far too much peaked to match the mor- 
phological shape of the observed one. Figure 2 illustrates that it shows a slightly better 
resemblance to the Holmberg et al. (2007) distribution for a spherical solar neighborhood 
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Figure 2. G-dwarf metallicity distribution with two-infall model for SD+DD (gray histogram) 
Spherical solar neighborhood observations by Holmberg et al. (2007) (black histogram). 
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as is used for comparison by Kobayashi & Nakasato (2011). However, since in our code 
a cylindrical model is used, we should indeed compare to such observations. 

The metallicity distribution which is found when a constant SFR (with total star 
formation equal to the previous case) is assumed, is shown in Fig. 3. In addition to cylin- 
drical model observations, this figure shows the histogram obtained with the combined 
SD (with ci = 1) + DD model, as well as the models for only SD (with c\ — 3) and DD 
respectively. It shows that the combined model produces a very good match, while either 
scenario alone (especially the SD, even with high ci) is unable to produce sufficiently high 
metallicity. Therefore, this study with its updated and internally PNS-computed SN la 
yields supports previous results favoring the SD + DD model. This combined model also 
satisfactorily reproduces the observed [C/Fe] and [O/Fe] vs. [Fe/H] relations. 

4. Conclusions 

We find that both the delay time (see also Mennekens et al. (2010)) and G-dwarf 
metallicity distribution point toward a significant contribution by both the single and 
double degenerate scenario, the latter mainly through systems that have undergone a 
quasi-conservative Roche lobe overflow followed by a common envelope phase. The critical 
dependence of both distributions on certain binary evolutionary processes, exactly those 
processes which are described in population codes by still uncertain parameters, is a way 
to find out more about these processes and thus further constrain these parameters. 
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Figure 3. G-dwarf metallicity distribution with constant SFR for SD+DD (gray histogram), as 
well as for SD (gray line) and DD (black line) only. Cylindrical solar neighborhood observations 
by Holmberg et al. (2007) (black histogram). 




